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Abstract. We study the singularity formation of smooth solutions of the relativistic 
Euler equations in (3 + l)-dimensional spacetime for both finite initial energy and 
infinite initial energy. For the finite initial energy case, we prove that any smooth 
solution, with compactly supported non-trivial initial data, blows up in finite time. 
For the case of infinite initial energy, we first prove the existence, uniqueness and 
stability of a smooth solution if the initial data is in the subluminal region away from 
the vacuum. By further assuming the initial data is a smooth compactly supported 
perturbation around a non-vacuum constant background, we prove the property of 
finite propagation speed of such a perturbation. The smooth solution is shown to 
blow up in finite time provided that the radial component of the initial “generalized” 
momentum is sufficiently large. 


1. Introduction. 

In this paper, we study the singularity formation of solutions of the Einstein 
equations for an isentropic perfect fluid. Due to the hyperbolic nature of these 
nonlinear equations, one expects singularity formation in the solutions. Indeed, 
one even expects black holes to form. However, singularity formation in relativistic 
flow is not yet well-understood; the theory is most lacking in the multi-dimensional 
case, (3 + l)-dimensional spacetime. 

As a first step in this direction, we consider here the relativistic Euler equations 
for a perfect fluid in 4-dimensional Minkowski spacetime, 

DivT = 0, (1.1) 

where 


T ij = (p + pc^vdu 3 + pg ij , (1.2) 

is the stress-energy tensor for a perfect fluid, and g 13 denotes the flat Minkowski 
metric, g 13 — diag(— 1 , 1 , 1 , 1), x = (x°, x 1 , x 2 , x 3 ) T with x° = ct. p is the mass- 
energy density, p is the pressure, c is the speed of light, and u is the 4-velocity of 
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the fluid. Recall that since u = 


Minkowski space), it follows that 


1 dx 
c dr 


(t is the proper time, u is a unit 4-vector in 


(« 0 ) 2 - X >“) 2 = i . 

a=l 

and thus only three of the quantities u°, u 1 , u 2 , u 3 are independent. We now fix our 
space-time coordinates as (t, x 1 , x 2 , x 3 ) T , set x — (x 1 ,x 2 ,x 3 ) T , u = (u 1 ,u 2 ,u 3 ) T , 
and let 

cu 

v= y(i + M 2 )' 

One easily derives from equation (1.1) the relativistic Euler equations: 


a , P C +P 
1 c 2 __ y 2 c 


P_ 

2 - 


+ V a 




c 2 — V 2 


dt{ ei±p v) + Vx 


(d^-^v®v)+V lP = 0, 


(1.3) 


in the unknowns p, v and p. Here V x denotes the spatial gradient operator. Given 
a scalar k and 3-vectors a and b, by the notion a®b we mean the matrix ab T , 
while 


V x • (kab T ) = (Va, • (fcaib), V x • (fca 2 b), V x • (fca 3 b)) T . 


We consider the Cauchy problem for (1.3) with initial data 


p{ x i 0) = po(x), v(x, 0) = vq(x). (1.4) 

Equations (1.3) close if we assume an equation of state, p — p(p), p{ 0) = 0 with 

P(P) > 0, 0 < p'(p) < c 2 , p"(p) > 0, for p E (p*, p*), (1.5) 

where 0 < p* < p* < oo. For a 7 -law, p(p) = cr 2 p 7 with 7 > 1, the constant p* is 
chosen as follows: if 7 = 1, then p* = 00 ; and if 7 > 1, then p'(p*) — c 2 . Thus the 
unknowns for the Cauchy problem (1.3)-(1.4) are p and v. For more details on the 
derivation of equations (1.3) and a discussion of (1.5), see [14]. 

We are interested in the life span of smooth solutions for the Cauchy problem 
(1.3)-(1.4). For this purpose, we shall discuss two different cases: the case of finite 
initial energy, and the case of infinite initial energy. For the first case, we shall prove 
that if the initial data has compact support, then the life span of any non-trivial 
smooth solution for the Cauchy problem (1.3) and (1.4) is finite. For the second 
case, we show that if the initial data is a compactly supported perturbation around 
a non-vacuum background, then the life span of smooth solutions is finite provided 
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that the radial component of the initial “generalized” momentum is sufficiently 
large; c.f. Theorem 3.2. 

We start with the infinite energy case. The local existence of classical solutions 
of the Cauchy problem (1.3)-(1.4) has been established by Makino and Ukai ([6, 7]) 
provided that the initial data is in the subluminal region away from the vacuum. 
A sharper result is proved here in Theorem 2.1 in Section 2, where the stability of 
the solution with respect to the initial data (c.f. Corollary 2.2) and the properties 
of finite propagation speed (c.f. Lemma 2.3) are presented. In Section 3, we first 
derive some interesting structural properties of (1.3) in Lemma 3.1, then we prove 
a blowup result for smooth solutions (c.f. Theorem 3.2). Our proof is in the 
spirit of the work of Sideris for classical Euler equations [11] and is based on the 
largeness of the initial radial component of “generalized” momentum, which of 
course implies the largeness of the initial velocity. However, in our case, the velocity 
is still subluminal. In Section 4, we prove our blowup result for smooth solutions of 
(1.3)-(1.4) with lion-trivial initial data that has compact support. In Section 5, we 
make some remarks concerning our results. A discussion on the type of singularity 
is also given. The existence of initial data satisfying our blowup conditions is also 
shown there. All the results in Section 2 are based on the existence of a strictly 
convex entropy function for (1.3), which was constructed by Makino and Ukai in 
[6, 7]. For the reader’s convenience, we present the construction in the Appendix, 
correcting a few errors in the original papers. 

Before proceeding, we now briefly review the methods and results of singularity 
formation for nonlinear hyperbolic systems. In one space dimension, the theory is 
fairly complete. It was proved that a singularity develops in finite time no matter 
how small and smooth the initial data is; c.f. [4, 5, 13]. These results were estab¬ 
lished by the characteristic method, which is quite powerful in one space dimension. 
In more than one space dimension, there are no general theorems available mainly 
because the characteristics become intractable. However, the approach via certain 
averaged quantities was introduced by Sideris [11] to prove the formation of singu¬ 
larities in three-dimensional compressible fluids. This idea avoids the local analysis 
of solutions. A similar technique was used to prove other formation of singularity 
theorems. We refer to [8], [9], [16] for classical fluids, and [2], [10] for relativistic 
fluids. Blowup results for relativistic Euler equations are announced in [2] and [10]. 
However, as remarked on page 154 of [2], “ the unpublished proof in [10] contained 
an error which invalidated the argument ”. Furthermore, we note that the coeffi¬ 
cient matrices in (2.15) of [2] constructed through (2.16) of [2] are not symmetric 
away from the equilibrium. But the symmetry of (2.15) in [2] is crucial to prove 
the finite propagation speed property needed in their proof. Thus the argument in 
[2] is not complete. Furthermore, we note that the equation of state used in [2] and 
[10] is different from ours. In addition, the approach of [2] is also different from 
ours. Our approach is closer to the method of Sideris [11], Finally, we remark that 
the equation of state (1.5) in this paper is interesting for cosmology. It includes 

3 



many physical cases, e.g. 7 -laws, p(p) = cr 2 p 7 , 7 > 1. For instance, the case 

p(p) = \ c2 P 

is very important in cosmology; it is the equation of state for the Universe in 
earliest times after the Big-Bang; see [15]. Some cases discussed in [2] (e.g. when 
s — const.) satisfy (1.5) as well. Another important example (see [15, p.319]) is the 
equation of state for neutron stars, where 


P = Ac 5 a{y ), p = Ac 3 b{y ), 

py 4 ry _ 

a (v) = J ^i + c 2 dq ’ b ^ = 3 J q 2 V 1 + q 2 <k- 


Here A is a positive constant. This equation of state implies the following asymp¬ 
totics: p —> \c 2 p as p —> cx) and p —> |A 2 / 3 p 5 / 3 as p —> 0. It is easy to see 
that 


P'(p) 


3(1 + y 2 ) 


> 0 , p"(p) 


2 

9 Acy 


(1 T y 2 )—^/ 2 > 0 , 


whenever, y > 0. We also note that y = 0 is equivalent to p = 0. Thus the equations 
(1.6) also satisfy (1.5). 


2. Existence of Solutions: Infinite Energy Case. 

In this section, we consider the local existence of smooth solutions for the Cauchy 
problem (1.3)-(1.4) when the initial data is away from the vacuum. For this pur¬ 
pose, we introduce some convenient notation: 


P = 


pc 2 + p 


. pc +p 
P ( — 2 2 " 

C — V A 

The Cauchy problem (1.3)-(1.4) becomes 



( 2 . 1 ) 


f Pt + V* • (pv) = 0 

< (pv) t + V x • (pv ® v) + V x p{p) = 0, (2.2) 

[ p(x, 0 ) = p 0 {x), v(x, 0 ) =v 0 (x). 

Let p* < p* be non-negative constants in (1.5) subject to the subluminal condi¬ 
tion p'(p*) < c 2 . We set 

2 = {p,v 1,V2,V 3 ) T 

and define the region il z by 

O z = {z : p* < p < p*, v 2 < c 2 }. (2.3) 
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Theorem 2.1. Assume an equation of state is given as in (1.5). Suppose the 
initial data zq(x) = (po(x), vq(x)) t is continuously differentiable on R 3 , taking 
values in any compact subset D of 0 2 and that V x zo(x) G H l ( R 3 ) for some l > 
3/2. Then there exists Tx, 0 < T x < oo, and a unique differentiable function 
z(x, t ) = (p(x, t), v(x, t)) T on R 3 x [0, T*,), taking values in FL Z , which is a classical 
solution of the Cauchy problem (l.S)-(l.f) on R 3 x [0,Too). Furthermore, 

V x z(;t)eC°([0,T oo )-H l ). (2.4) 

The interval [0, Tx,) is maximal, in the sense that whenever Tx, < oo, 

lim sup || V x z(-, = °° (2-5) 

t > T 00 

and/or the range of z(-,t) escapes from every compact subset of Fl z as t —> Tx,. 

This theorem will be proved by applying Theorem 5.1.1 in Dafermos [1] for 
hyperbolic conservation laws endowed with a strictly convex entropy. We state this 
theorem here for readers convenience. 

Theorem A. Assume that the system of conservation laws 

m 

Ut + Y. d *« G *(U) — 0, x E R m , T G O c R n , (*) 

(X=l 

is endowed with an entropy r) with W 2 rj{U) positive definite, uniformly on a compact 
subset of O. Suppose the initial data U(x, 0) = Uq(x ) is continuously differentiable 
on R m , takes values in some compact subset of O and VUq G Ft 1 for some l > m/2. 
Then there exists T 00 , 0 < Tx < oo, and a unique continuously differentiable 
function U on R ?n x [0, Tx), taking values in O, which is a classical solution of the 
initial-value problem (*) with initial data Uq on [ 0 ,Tx). Furthermore, 

\/Uf,t)eC°([0,T oo ) ] H l ). 

The interval [ 0 , Tx) is maximal, in the sense that whenever Tx < °° 

lim sup || VC/ (•, t)|| T°o — OO 

t 

and/or the range of U(-, t) escapes from every compact subset of O as t —> Tx- 

Proof of Theorem 2.1. We first rewrite (1.3) or (2.2) in the form of conservation 
laws, 


», + £(/*("))» = 


k =1 


( 2 . 6 ) 
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(2.7) 


where 6 = (0 O , 6 U 0 2 , 0 3 ) T and f k (9) = (9 k , /f, / 2 fc , / 3 fc ) T are defined by 


= P: Oj = pVj, 

fj = P v f ,J k + J = 1, 2, 3. 


By Theorem A, it is sufficient to show that (2.6) has an entropy r](9) with X7 2 p(9) 
positive definite in Vt z . Such an entropy, due to Makino and Ukai [7], is constructed 
in the Appendix of this paper. 

Define 


<Kp) = 


' Pm 


rc 2 + p(r) 


dr, K — p m c 2 +p(p m ), 


( 2 . 8 ) 


p m being any fixed number in (p*, p*). The entropy given in (6.25) below is 


V 


c 2 p- 


cKe^^ 
\J c 2 — v 2 


(2.9) 


We now verify that V 2 p(0) is positive definite in Vt z . To this end, we first compute 
Vep(9). By the chain rule, we have 


w T = (V d y) = (,V z ri)(V z e)-\ 


where (V z 9) 1 is defined in (6.10), and w T = (wo,Wi,W 2 ,w 3 ) is given by 


Wn - __cW(p) , 2 

U 0 ^ c 2_.y2p/2 T C > 

_ c<K(p) 


W 


, . — . n — 1 2 *1 

] ( C 2_ v 2y/2 U ]1 J 


( 2 . 10 ) 


with 


d>(p) 


(pc 2 + p ) 


( 2 . 11 ) 


We remark that w can serve as a symmetric variable which reduces (1.3) to a 
symmetric hyperbolic system [1, 3]. For the Hessian matrix H of rj, we compute 


H =V 2 p(60 = V e w T = (V z w T )(XZ z 9)~ l 

= c$(p)£i 

(pc 2 +p)(c 2 -V 2 ) 1 / 2 1 

= c$(p)^! / A! A 2 n T \ 

(pc 2 Tp)(c 2 — v 2 ) 1 / 2 \ A 2 v A 3 vv T + A 4 I 3 y 

Here, is given in (6.11) below, and the A t are given by 


A\ = c 4 (p'c 2 + 2 p'v 2 + c 2 v 2 ), A 2 = — c 2 (c 4 + 2 c 2 p + p'v 2 ), 

As = (c 4 + 2 c 2 p' + p'v 2 + 2 p'(c 2 — v 2 )), A 4 — (c 2 — v 2 )(c 4 — p'v 2 ). 
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( 2 . 12 ) 


(2.13) 



We now show that H is positive definite. From (2.12), we see that it is sufficient 
to show H i is positive definite. Let r = (ro, r T ) T be any 4-vector with r G R 3 . We 
calculate: 


H t r = (r 0 , r ) 


( Ai A 2 v T 
l A 2 v A 3 vv t + A4I3 


(r 0 ,r) 


= (A t rl + 2 A 2 r 0 v T r + A 3 (v T r) 2 + H 4 r 2 ). 


Letting Ai — (1 — 5)Hi with | > 5 > 0 to be determined in (2.14) below, we have 


(AiTq + 2A 2 r 0 v 1 r + A 3 (v T r) 2 + A A r 2 ) 


A 2 t \ 2 / 1 / 10 , , \ 5 A3) 


= Ai(ro + ^'rf - {—{A 2 - A X M) + . * A 

A 1 A t 1 - d At 


2 ){y T r ) 2 + SAtr'y + A^r z 


1 


2 $ ^ 2 „. 2\„2 


> (A t - j-(A\ - A,A 3 )v l _ 6M 


v 2 )r 2 + 5AtrQ 


> (Lw; 


p'(c 2 -V 2 ) 2 (c 4 -p'v 2 ) Al 2,2 


(p'c 2 + 2 p'v 2 + c 2 v 2 ) At 

> SAtrQ + 5r 2 . 


2S-Av 2 )r 2 + SAtr 2 


Here, we determine <5 by 


0 < S + 25^v 2 < 
At 


p'{c 2 — u 2 ) 2 (c 4 — p'v 2 ) 
(p'c 2 + 2 p'v 2 + c 2 v 2 ) 


We thus conclude that 


(2.14) 


r T iL!r > (SAy + 5r 2 ). 

This proves H 1 is positive definite in Q~. Hence, H is positive definite, and r/ is 
strictly convex on Vt z . This completes the proof of Theorem 2.1. □ 

The existence of a strictly convex entropy guarantees that classical solutions of 
the initial-value problem depend continuously on the initial data, even within the 
broader class of admissible bounded weak solutions; see [1]. Here, by admissible 
bounded weak solution, we mean bounded functions satisfying the initial value 
problem and entropy inequality in the sense of distributions. The following Theorem 
B is Theorem 5.2.1 in Dafermos [1]: 

Theorem B. Assume that the system of conservation laws (*) is endowed with 
an entropy p with V 2 p(U) positive definite, uniformly on compact subset of O. 
Suppose U is a classical solution of (*) on [0,T), taking values in a convex compact 
subset N of O, with initial data Uq. Let U be any admissible weak solution of (*) 
on [0, T), taking values in N, with initial data Uq. Then 
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I \U(x,t) — U(x,t)\ 2 dx < ae bt / \Uq(x) — Uq(x)\ 2 dx 

J\x\<R J\x\<R+st 

holds for any R > 0 and t E [0, T), with positive constants s, a, depending only on 
N, and a constant b that also depends on the Lipschitz constant ofU. In particular, 
U is the unique admissible weak solution of (*) with initial data Uq(x) and values 
in N. 

The following corollary is a consequence of the convexity of rj and Theorem B. 

Corollary 2.2. Let 6 be a classical solution of (2.6) obtained in Theorem 2.1 with 
initial data 9o(x ) taking values in a compact subset D of Q z , and let 0 be any 
admissible weak solution of (2.6) on [0, Too), taking values in D , with initial value 
9q(x) E D. Then 


\9(x, t ) — 9(x, t)| 2 dx < ae 


bt 


\9q(x) — 9q(x)\ 2 dx 


|ie|<-R 


I as I <R+st 


holds for any R > 0 and t E [0, Tx,), with positive constants s, a, depending only on 
D, and a constant b that also depends on the Lipschitz constant of 9. In particular, 
9 is the unique admissible weak solution of (2.6) with initial data 9 0 (x) and values 

in D. 

In the next lemma, we will show that a compactly supported perturbation around 
a non-vacuum background propagates with finite speed. For this purpose, we con¬ 
sider the following Cauchy problem 


Pt + • (pv) = 0 

{pv) t + V x • (pv ( 8 ) v) + V x p(p) = 0, 

(po(x) - p, v 0 (x)) = 0, for \x\ > R, 
where, R > 0, and 0 < p E (p*, p*) satisfies the subluminal condition, 


(2.15) 


p'(p*)<c 2 . (2.16) 

Lemma 2.3. Let (p,v)(x,t) be a C 1 solution of the Cauchy problem (2.15) (equiv¬ 
alent to (1.3) with the same initial data), where the initial data ( Pq,vq)(x ) takes 
values in a compact subset ofVt z (c.f. (2.3)). Then the support of (p — p,v)(x,t) is 
contained in the ball B(t) — {x : |x| < R + st} where 

s = y/p'(p) 


is the sound speed in the far field. 

Proof. This lemma is a consequence of the local energy estimates. It will be proved 
using the method of [11] for symmetric hyperbolic systems. For this purpose, we 



first observe that w = (V gp) T given in (2.10) (where r/ is as in (2.9)) renders the 
system (2.15) symmetric hyperbolic [3]: 


A^ + ±A^ = 0, ,2,7) 

i=1 

where the coefficient matrices A a (w ) = (A^ n ), a,m,n = 0 , 1,2,3 are given by 
(6.5); that is 

A° = (v^r\ ^(v^xvfor 1 . 

We now compute the explicit form of these matrices. First of all, we have 


= (Vjv)- 1 = 9(p) 


a i a 2 f 

a 2 f a 3 vv T + 04/3 


where /3 is the 3x3 identity matrix, and 


^(p) =-^(pc 2 +p) 2 e 


and 


/„,2 


ai = 


«3 — 


c + 3p n 


c 4 + 2p'c 2 + p'v 2 
c 3 p’{c 2 — v 2 ) 3 / 2 

1 

c(c 2 — U 2 ) 1 / 2 


c 3 p'(c 2 — V 2 ) 3 / 2 


1 a 2 — 


c 2 + 3p ; 


cp'(c 2 — v 2 ) 3 / 2 
For A fc , k = 1, 2, 3, we first compute 

(V,/ fc ) - (VJ^VJ )- 1 

0 


, 04 = 


[c 2 (c 2 + v 2 )p'Eie k - c 2 (c 2 + p^ExVkv] [-C 4 ve£ + v k I 3 ] 
with e fc = (hifc, 5 2 fc, <W T , where (VJ fc ) is given by 

(V z / fc ) = 

c.f. (6.13) below. Then we have 


B 3 v k B 2 v k v T + B 4 e k 

B 3 v k v + p'e k B 2 v k vv T + B 4 ve k + B 4 v k I 3 ) ’ 


(2.18) 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


A k 


(V */*)^)" 1 = (V e f k )A° 

\ f °2 V k [a 3 v k v T + a 4 el\ 

' ' V [ fl 3 v kV + a 4 e k ] [a 3 v k vv T + a 4 v k I 3 + a 4 (e k v T + ve[)] 
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( 2 . 22 ) 



It is clear that the matrices A k (w ) are all real symmetric and smooth in Q,. Fur¬ 
thermore A Q {w) — (V 2 ??) -1 is positive definite in 

Now, we choose p m = p subject to (2.16) for convenience. In this setting, the 
background state in the re-variable becomes ui = w(p, 0) = 0. Set 

A\w ) = (A°{w))~ 1 A i {w), 


and define 


3 

Q( A,O = A/4-^?«i“(0). 

C£=l 


(2.23) 


where (A, ()eRx5 2 ( S 2 is the unit 2 -sphere). Using the real symmetry of A a ( 0), 
for each £ G S' 2 , we see that the characteristic equation 


detQ( A, £) = 0, 

has real roots A*(£), i = 0 ,1,2,3, called the characteristic speeds. Let A be the 
largest absolute value of these characteristic speeds. For any fixed (xo,to) G R 3 x 
( 0 , Too), we define the family of cones 

C T — {(s, x) : |x 0 — x\ < X(to — s), 0 < s < r}, (2.24) 

parametrized by r G [ 0 , to) and the associated cross sections 

E a = {{fi,x) G C T : n = a}. (2.25) 

We introduce the linear partial differential operator 

P = A »W t +t Aa (0)A r , (2.26) 

a=l 


where again 0 is the background state in the re-variable. The equation (2.17) reads 


Pw = (n°(0) - A°(w))— + Y,(A a (0) - A“(w))- 


ot=l 


dXr 


Now we multiply both sides of (2.27) by 2W 1 , to get 

3 8 

d t [w T A 0 (0)w\ + ^ ~— [w t A°‘(0)w\ 

= 2w t 


Ct = l 


(A°(0) - A°(w))?£ + 5>«(0) - A°(w)) OW 


ot =1 
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dx r 


(2.27) 


(2.28) 



We integrate (2.28) over C T to obtain 


[. 3 

/ d a [w T A a (0)w\(x, a) dxcla 

J Cr a = 0 

r 3 

= / 2«; T ^(A“(0)-W>))d> Q 


w 


(2.29) 


a=0 


< (7max|Vre| 

C T 


\w\ 2 (x, a) dxda. 


'0 JE a 


Where, do — dt, and we have used the mean value theorem in the estimation of the 
last step. For the left hand side of (2.29), we want to apply the divergence theorem 
since it is in divergence form. For this purpose, we need to determine the boundary 
of C T and the associated unit outer normal vector. The boundary of C T consists of 
three parts: the cap E r with unit outer normal (1,0, 0,0) T , the base Eq with unit 
outer normal (—1, 0, 0, 0) T , and along the surface 


the unit outer normal vector is 

1 

n — 


i — x\ — X(to — cr),0 < a < r} 

(2.30) 

(X J .T\T - (x 0 -x) 

A ’ ¥ ’ \x-xA ' 

(2.31) 


\/l + X 2 

To see (2.31), we note that on R r one has 

(A, -v T ) T • ((t 0 - cr), (x 0 - x) T ) T = 0. 

We now apply the divergence theorem to the left hand side of (2.29), 

3 

/ d a [w T A a (0)w](x, t) dxdr 

J Cr a = 0 


(re T A 0 (0)«;)(a:, r) dx 


(w t A°( 0)w)(x,0) dx 


'e t 


'E 0 


(2.32) 


+ 


Va 2 + i 


'0 JdE. 


(Xw t A°(0)w — w T u a A a (0)w)(x, a) dS x da, 




where dS x denotes the surface element on dE a . The third term of (2.32) on the 
right hand side can be simplified as follows: 

3 

(Xw T A°(0)w — w T ^ E a A°‘(0 )w) 

a=l 


3 

= w T A°( 0 )(X /4 — v ot A a (0))w 

ct =1 

= m; t A 0 (0)(5(A, v)w. 
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(2.33) 



We recall that A°(0) > 0 and 


3 

A°(0)Q(A, v) = AA°(0) - WA1“(0) (2.34) 

a=l 

is real symmetric. We claim that for any v G ,S' 2 . 

A°(0)Q(A, v) > 0, (2.35) 

which will be verified at the end of this proof. 

Therefore, we conclude from (2.29), (2.32), (2.33) and (2.35) that 


(w T A°(0 )w)(x,t) dx 


'Et 


(2.36) 


< / (wr A u (0)rc)(a;, 0) dx + Cimax| Ww\ 

Je 0 w 


\w\ 2 (x, a) dxcla. 


'0 JE a 


Since A°(0) > 0, there are positive constants and C$ such that 


/ \w\ 2 (x,r) dx < C2 / |u>| 2 (x, 0) dx + C3 / / \w\ 2 (x,a) dxda, 

J E t J Eq Jo jEcr 

which, by Gronwall’s inequality implies that 


/ \w\ 2 (x,r) dx < C f 2 e c ' 3r ( / |rc| J (a;,0) dx). 

' E t Je 0 


(2.37) 


Therefore, if w(x, 0) = 0 for \x — Xq\ < Xto, then w(x, r) — 0 for any r G [0, to) and 
\x ~ Xq\ < A (to — r). This implies that, if w(x, 0) = 0 for |x| > R, then w(x, t) = 0 
for |x| > R + Xt. 

The next step is to verify that A = y/p'(p). For this purpose, we compute 
the largest possible characteristic speed at a constant background state. Now we 
compute the eigenvalues of 


3 

]T6*i Q (0), £gS 2 . (2.38) 

a=l 

Since 

A a = (V 2 r ? )(V,D(V 2 r ? )- 1 , 

the matrix in (2.38) is similar to the matrix 

3 

M(() = Y, W, ? 6 S 2 , (2.39) 

a=l 
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where 6 = (p, 0,0,0) 7 is the background state in the 6- variable. It is easy to 
compute: 

v e f a (e) = v z f a (v z ey 1 \ p =- P!V =o 

= ( 0 e n (2.40) 

\ v p'(p)e« 0 ) ’ 


where e* = (hii, , $ 3 i) T ■ Thus, one has 


M(0 




(2.41) 


Now we claim that M (£) has rank 2. This can be seen by 


MM t 


(o p'(p) 2 ^ T )’ 


since rank(^ T ) — 1. Thus, M(£) has two non-zero eigenvalues and two zero 
eigenvalues. We need to find all non-zero eigenvalues of M(£). We compute 


0 = det(M(0 - r/ 4 ) 


det f ** 

VP (p)£ - r/ 3 

dCt ((p'(p) -^ 2 )^ 

(p'(p) -r 2 )det Q 



(2.42) 


and this implies ±a/ p'(p) are the two distinct non-zero eigenvalues of M(£). There¬ 
fore, we have 

A = V7W). (2.43) 

Notice that we did not use (2.35) to obtain (2.43). 

The last step is to verify (2.35). Since A = s — yV(p), we have 


A » = l(‘o I ) • ^“(0) = 4 ( 0 e “l,a= 1,2,3, 


c 2 \e a 0 


and thus 


a 0 (0)Q(a,O) = 4 

C 2 
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-// sT 


(2.44) 



Let r = (ro, r T ) J be any 4-vector with r G R 3 , we compute: 


r T H°(0)Q(A, u ))r = — (r 0 , r T ) 


-v sh 


(ro, r 1 


Here, we have used 


= -w(- 7 q - 2 r 0 v r + slrl ) 
c z s 


> -w(-rg - 2 r 0 |r| + s|r| ) 

C z S 


= —y(r 5 - 2 sr 0 |r| + s 2 |r| 2 ) 
sc z 

= (r o - s|r |) 2 > 0 . 


( is T r ) < \J (z/ 2 )(r 2 ) = |r|. 


(2.45) 


Thus, we have proved (2.35). The proof of this lemma is complete. □ 

3. Singularity Formation: Infinite Energy Case. 

In this section, we prove the singularity formation of smooth solutions of (2.15) 
when the initial radial “generalized” momentum is large. To begin, we prove the 
following two easy but useful identities. 

Lemma 3.1. p and p satisfy the following identities: 

, 1 _ 2 

P = ~P V +P> 

i „ < 3i > 

p=—pv + <p+ 
c c 

Proof. From (2.1), it is easy to see 


. rpC Z +P p 

p = “2-2 “ "2 

c z — v z c z 

pc 4 + pc 2 — pc 2 + pv 2 
c 2 (c 2 - V 2 ) 
pc 2 c 2 + pv 2 
c 2 (c 2 - V 2 ) 

pc 2 v 2 + pc 2 (c 2 — v 2 ) +pv 2 
c 2 (c 2 -v 2 ) 

pc 2 + p 2 pc 2 (c 2 —v 2 ) 

c 2 (c 2 -v 2 ) V + c 2 (c 2 — V 2 ) 
1 ~ 2 

= — pv + p. 
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Hence, p = p + implies 

1 _ o P 

p = —pv + p+—. 
c z c z 

We again denote the sound speed in the far field by 

« = vV(p)> 

and define the following quantities: 

M(t) = j [p(p, v) - p(p, 0)] (x, t ) rix, 
Fit) = jpv. X d X . 


□ 


(3.2) 


By Lemma 2.3, both M{t) and F(t) are well-defined as long as the smooth solution 
exists. Using these two quantities, we shall show that the smooth solution of (2.15) 
obtained in Theorem 2.1 blows up in finite time if the initial data is subject to some 
restrictions. Roughly speaking, if M(0) > 0, and -F(O) > 0 is sufficiently large, then 
the solution will blow up in finite time. 

Theorem 3.2. Assume that the initial data of (2.15) and p are chosen such that 
M(0) > 0, F(0) > 0 and s 2 < ±c 2 . If 


F{ 0) > T 


327 rs 



R 4 maxpo(ic), 


then the smooth solution of the Cauchy problem (2.15) obtained in Theorem 2.1 
blows up in finite time. 

Proof. Using Lemma 3.1, we know that 


thus, 

This implies that 


P(p,v) = —pv 2 +p, 
c z 

p(p, 0) = p. 


p-p = (p(pw) - p(p,0)) - 
15 



(3.3) 

(3.4) 



From the first equation of (3.2), it is easy to see that 


d f 

^M(t) = J (p{p,v) - p(p,0)) t dx 
= [ p t dx 


V x • (pv) dx 


= 0, 


where we have used the first equation in (2.2) or (2.15). Hence, 

M(t) = M{ 0) = j ( po(x ) — p) dx > 0. 
Using the second equation in (2.2) and (3.2), we compute 
F'(t) = [ ( pv)t • x dx 


= — J [V • (pv (8) v) + Vp(p)] • x dx. 

But if p — p(p), we have 

Vp • x — V(p — p) • x 

— V • [x(p — p)\ — (p — p) V • x 
— V • [x(p — p)\ — 3 (p — p). 

We also note that 


(3.5) 


(3.6) 


(3.7) 


3 

V • (pv ® v) • x = ^2 d Xi (pViVj)xj 
i,j =1 
3 

= ^2 (P v i v j x j ) - PViVjd Xi Xj\ 
i,j =1 

3 

= ^2 [dxiipviVjXj) ~ pviVjdij] 
i,j =1 

3 

= —pv 2 + ~22 V • (pu (E> vxj), 

3 = 1 


(3.8) 


where u 2 = u T u. Inserting (3.7) and (3.8) into (3.6), and using the divergence 
theorem, we obtain 
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F'(t) = / pv 2 dx + 3 {p — p) dx. 


( 3 . 9 ) 


Since p"(p) > 0, p'(p) is a non-decreasing function of p. It is clear that 


p(p)-p{p)= p'(0 >p'(p)(p~p)- 


(3.10) 


Thus, using (3.4) and (3.5), one has 


F'(t) > / pv 2 dx + 3s 2 / (p — p) dx 


3s 2 

= / pv 2 dx + 3 s 2 M(t) -y / pv 2 dx 


= ( 1 - 


3s" 


pv 2 dx + 3s 2 M(0) 


> (1 - ^-) y '/ 91 ’ 2 d x - 


On the other hand, we have the following estimate: 


F 2 {t) — ( pv • x dx)‘ 


(3.11) 


' B{t ) 

where we have used the following fact: 


— ( / |x| 2 p dx){ / pv 2 dx) 

J B(t) J B(t) 

< 2( / \x\ 2 p dx)( / pv 2 dx), 

J B(t) J Bit) 


(3.12) 


P < 2 P- 


(3.13) 


To see this, we note that from the subluminal condition p'(p) < c 2 , together with 
p( 0) = 0, we get p(p) < c 2 p. Thus 

p = \pv 2 + p + < \pv 2 + 2/9 = p + p < 2/3. 

c c c 


Due to (3.3) and (3.5), we have the following estimate 
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/ p\x\ 2 dx < (R + st) 2 / p dx 

' B(t) J B(t) 

= (R + st) 2 (M(t) + / p dx) 
J B(t) 

= (R + st) 2 (M( 0) + / p dx) 

J B(t) 

— (R + st) 2 / Po( x ) dx 
J B(t) 

47T c 

< — (R + st) (maxp 0 (x)). 
o 


Hence, (3.12) gives 


where 


F 2 (t) < — (R +st) 5 (max p 0 (x))( f pv 2 dx) 

3 JB(t) 

= K 0 (R + st) 5 ( / pv 2 dx), 

J B(t) 


K 0 = ^-(maxpo(x)). 


Thus (3.11) and (3.15) imply that 


3s 


nt)> (i-^)K oHr+ st)- 6 ^), 


„2 ' 0 


SO 


> KxiR + st)- 5 , 


F 2 


where Ki = (1 — ^-)/i 0 1 . Integrating (3.18) with respect to t, one has 


< —— ~ ^r-[R 4 - (R + st) 4 ] = ip(t). 


F(t) ~ F( 0) 4s 


Now ip(0) — > 0 by assumption, and 


if 


^ + °°> =m - % R " < “• 


4sT? 4 

F(o) > = r. 
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(3.14) 


(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 



Therefore, 

—= 0, for some t 0 > 0. (3.21) 

r (to) 

Thus the life-span T of smooth solutions satisfies T < to- This completes the proof 
of Theorem 3.2. 

□ 

4. Singularity Formation: Finite Energy Case. 

Due to the hyperbolic nature of Einstein equations, one expects the finite prop¬ 
agation speed of waves in the solutions. We will prove in the following lemma that 
for any smooth solution with compactly supported initial data, the support of the 
solution is invariant in time. 

Lemma 4.1. Let (p,v)(x,t) be a smooth solution of the Cauchy problem (1.3)- 
(1-4) U P to some time T > 0. If the support of initial data is contained in the 
ball Br( 0) centered at the origin with radius R, then the support of (p,v)(x,t) is 
contained in the same ball B r(0) for any t G [0, T). 

Proof. Assume that the initial support of the solution is contained in a ball Br( 0), 
the support of the smooth solution will remain compact by the hyperbolic nature 
of the system (1.3). We denote by x(t;x o) the particle path starting at Xq when 
t — 0, i.e., 

d 

—x(t;xo)=v(x(t-,x 0 ),t), x(t = 0;x o ) = x o , (4.1) 

and by S p (t) the closed region that is the image of Br( 0) under the flow map (4.1). 
Hence, the support of the smooth solution of (1.3)-(1.4) will remain inside S p (t). 
Thus, fixing any xq on the boundary of Br( 0), we have Po(xq) = 0 and uo(:ro) = 0, 
and x(t;xo) is on the boundary of S p (t). Furthermore, 

^x(t; x 0 ) = v(x(t; x 0 ), t ) = 0, (4.2) 

due to continuity of v(x,t) and the fact that x(t;x o) sits at the boundary of the 
support of the solution. Therefore, x(t;x o) = Xq for any t G [0, T) whenever 
|xo| = R. Hence, S p (t) = Br( 0). This proves this lemma. 

□ 

Based on Lemma 4.1, we shall prove the following blowup result. 

Theorem 4.2. Suppose the support of the smooth functions (po(x), vq(x)) is non¬ 
empty and contained in a ball Br( 0) centered at the origin with radius R. Then the 
smooth solution of (1.3)-(l-4) with the initial data (po(x),vo(x)) blows up in finite 
time. 
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Proof. We first introduce the following functions, 



p\x\ 2 dx, F(t) 


pv • x dx, E(t) 


p dx. (4.3) 


Here, H(t) is the second moment of p, F(t ) is the total radial “generalized” momen¬ 
tum, and Eft) is the total “generalized” energy. These functions are well defined in 
the domain where the smooth solutions exist. Interesting relations between them 
can be obtained by the following calculations. 

Eft) is conserved, because using (2.2) one has 


We thus have 



V x • (pv) dx = 0. 


m = s(o) 


j po(x) dx > 0, 


for lion-trivial initial data. 
For H(t), we have 


(4.4) 


H'(t) = ^ j p t \x\ 2 dx 

= " /[V x • (pv)]\x\ 2 dx 

— pv • x dx 
= F(t), 

where we have used the relation 

[Va; • (pv)]\xf = V x • (pv\x\ 2 ) — pv • 2x. 

From the second equation in (2.2) and integrating by parts, we have 


H"(t) = F'(t) = / (pv)t • x dx 


= — J ([V* • (pv (8) v)] • X + (V xP • x)) dx 

= ( pv 2 dx + J 3p dx 

= C 2 ( (\pv 2 + -^) dx. 

Jb r { 0) c c 
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(4.5) 


(4.6) 



By Jensen’s inequality, we have 


'b r { 0) 


M dx = ( h fl 3)W(^ 

V 3 ; (^i? 3 ) 

^ i 4?r P 3 )v(^ B r(0)P dX ) 

> ( t R )p( {lfm) ), 


so (4.6) and (4.7) imply 


#"(£) > c 2 


[ \pv 2 dx + 4(-!r R 3 )p(pB ) 
Js R (o) c 2 c 2 3 


= c 2 iV(t), 


where 


_ Ib r ( o) Z 9 dx 
PB ~ (^-R 3 ) ’ 


is the mean density over 5 r( 0). Since 

i7(t) = i7(0) = / (4/w 2 + p) cte, 

JBr(O) c 

it is possible to bound iV(t) from below using E( 0). We consider two cases. 


'Bfl(O) C 


^ pv 2 dx > \ ) E(0 ), 

Zi 


we get 

iV(£)> 

On the other hand, if 

f 1 ~ 2 




•FBk(O) c 

then as 



1 

2 


£'(£) = E(0) = / (—pv 2 + p) dx < —E( 0) + / p dx, 

Jb r ( o) c 2 JBh(O) 


we have 


'Bh(O) 


p dx > —E( 0). 
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(4.7) 


(4.8) 


First if 


(4.9) 



Thus 


JV(t) > 4(hj? 3 )p(p B ) 

> ^(Y R3)pi l EBm 

(4.10) 

= SiU(O) > 0, 


where Bi = f i^ s y } p{\E B ({))), and E B { 0) = (J'Jsy Define B 
then (4.8)-(4.10) imply that 

= c 2 min{i,Bi}; 

H"(t) > BE( 0) > 0. 

(4.11) 

This gives a lower bound on H(t ): 


H{t) > -BE(0)t 2 + F(0)t + H( 0). 

2 

(4.12) 


In order to refine (4.12), we estimate F(t) in terms of H(t) and E(t). Using 
(3.10), we have 




Thus (4.12) or (4.14) together with (4.15) imply that the life-span of the smooth 
solutions must be finite if E(0) > 0. This completes the proof of Theorem 4.2. 

□ 

Notice that, for noil-trivial initial data, we have 


H( 0) 

_ 1 
” 2 
_ 1 
” 2 
< 0 . 


- ^E(0)R 2 

I po\x\ 2 dx — -R 2 I po dx 


'Bfl(O) 


Po{\x \ 2 — R 2 ) dx 


(4.16) 


This enables us to estimate the life-span as follows: from (4.14) and (4.15) we have 
for smooth solutions 

^BE{0)t 2 - D\JH(0)E(0)t + H( 0) < ]-E{0)R 2 . 

Z Z 

This is equivalent to 


(fit) = Bt. 2 — 2 Ddt + 2d 2 — R 2 < 0, 


(4.17) 


where d 2 = and 0(0) < 0 by (4.16). Hence, the life-span T of the smooth 

solution satisfies 

Dd + V D 2 d 2 - 2 Bd 2 + R 2 B 


T < 


B 


(4.18) 


5. Concluding Remarks. 

We have proved the blowup of smooth solutions of relativistic Euler equations in 
both cases: finite initial energy (Theorem 4.2) and infinite initial energy (Theorem 
3.2). In contrast to the characteristic method, we adapted the approach via some 
functions: total “generalized” energy, total radial “generalized” momentum, and 
the second moment. Our approach depends on the beautiful structure of the equa¬ 
tions and several quantities constructed from the natural variables; c.f. Lemma 
3.1. Although the relativistic Euler equations are much more complicated than the 
classical Euler equations, these structures make our proofs possible. We will now 
make some remarks on our results and discuss some related issues. 


Remark 1. In our blowup theorems, the velocity in the far field is assumed to 
be zero initially. For the more general case, say vq(x) — v off a bounded set, the 
change of variables (Sideris [11],) 


v —> v — V, x —> X + tv 

will reduce this problem to the case we considered. 
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Remark 2. The condition 

P\P) < j 

in Theorem 3.2 arises naturally in the proof. Here, 3 is the spatial dimension. In 
d dimensions, 3 is replaced by d. In particular, for d — 1, this condition is that the 
sound speed is subluminal. For p(p) — a 2 p, d = 1, this condition guarantees the 
genuinely nonlinearity of the relativistic Euler equations and allows the existence of 
global solutions in BV; see Smoller and Temple [14], This condition is not required 
in Theorem f.2. 

Remark 3. Our blowup results crucially depend on the compact support of the 
perturbations. Singularity formation for more general initial data remains open. 

Remark 4. The type of singularity which occurs is another open problem. The 
possibilities are: a) shock formation, b) violation of the subluminal conditions; e.g. 
|t>| tends to c, orp'(p) —> c 2 , c) concentration of the mass. 

For p(p) = a 2 p and d = 1, the singularity must be a shock if the initial data is 
away from the vacuum. It was shown in Smoller and Temple [14] that weak solutions 
exist globally in time with bounded total variation, subluminal velocity and positive 
density uniformly bounded from above and below. Furthermore, Smoller and Temple 
proved in [14] that the subluminal condition guarantees the genuine nonlinearity of 
the equations, so one concludes that the singularities in the solutions must be shocks 
by Lax’ theory [4], It would be interesting to clarify the types of singularities for 
relativistic Euler equations in multi-dimensions. However, black hole formation 
is impossible for our problem, since our spacetime is fixed to be flat Minkowski 
spacetime. 

Remark 5. The singularity in our Theorem 3.2 looks like shock formation. The 
largeness condition in radial “generalized” momentum, (3.20), implies that the par¬ 
ticle velocity must be supersonic in some region relative to the sound speed at in¬ 
finity. One can guess that the singularity formation is detected as the disturbance 
overtaking the wave front thereby forcing the front to propagate with supersonic 
speed. To see these things, we argue as follows. Using the fact p < 2 p(c.f. (3.13)), 
one has 


F( 0) < —i? 4 (max/5o(cc)) max |i>o(x)|, 
3 


(5.1) 


while 


4s 8f 4 

T = _ 3g2 y R maxpoO) 

C 2 

8n ^ , * / \ \ 

> 4s—it (maxpo(^)j- 

o 


(5.2) 


Hence, F(0) > T implies that 
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max |i?o (x) | > 4s. (5-3) 

This insures the initial particle velocity is supersonic in some region. However, the 
rigorous proof of shock formation is still open. 

Remark 6. The lower bound of the initial radial “generalized” momentum in 
(3.20) depends on the initial velocity through p. This is different from the Newto¬ 
nian case, where p is replaced by p and so in the Newtonian case it does not depend 
on the velocity. On the other hand, the velocity has to be subluminal. Therefore, we 
must show that the set of initial data required in Theorem 3.2 is non-empty. From 
(5.1)-(5.3), we find the set is non-empty if 

4<s 

c > max Iv 0 (a)I > --(5.4) 

A simple calculation shows that the necessary condition for s to satisfy is 

* < Cj- - |)c. (5.5) 

Since p > p, F( 0) is of the same order as the upper bound in (5.1). Thus, if s is 
chosen to be small (this can be done by choosing p small), one can easily find initial 
data satisfying the conditions required in Theorem 3.2. 

Remark 7. The equation of state p(p) satisfying (1.5) is quite general for isen- 
tropic fluids. It can be weakened by replacing p" (p) > 0 withp'(p) is non-decreasing. 
This includes the well-known 7 - law, p(p) = a 2 p 1 , 7 > 1 as a particular case. In 
fact, in the case of a 7 - law, (3.13) can be refined, and thus (3.20) can be replaced 
by a weaker condition, as we now show. 

When 7 = l, s = a. (3.13) is refined as p < (1 + ^r)p — p + ^ppv 2 by Lemma 
3.1. Thus, (3.20) can be weakened to: 

F{ 0) > Ti = 4<T 3a2 (1 + ^)^R 4 maxp 0 (x). (5.6) 

1 — —9— c o 

When 7 > 1, we observe that the subluminal condition 

p (p) — 7 a 2 p' y ~ 1 < c 2 

implies that ffi < ^p. Thus, 

„ R* , 1 1 . 2 

(1 + -)p = P + -P + —pv 

7 7 7 cr 

> p + -p 
7 

^ P 

>P+ — =P- 
c z 
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We can thus refine (3.13) to p < (1 + -)p, and then (3.20) is replaced by the 
following weaker condition: 


F{ 0) > r 2 


4- 1 4-7T ,, 

- -3^-a + -)^r R max po(x). 

1 _ 7 3 


(5.7) 


6. Appendix. 

For reader’s convenience, we justify the construction of a strictly convex entropy 
function for (1.3) due to Makino and Ukai in [7], and we will also correct several 
errors. To this end, we first record (2.6)-(2.7) here , 

3 

*< + £(/*(«))*„= 0 , ( 6 . 1 ) 

k=1 

where 9 = (0 O , 0i, 0 2 , # 3 ) T and f k (9) = (0 k . f(. f$, / 3 fc ) are defined by 

9o = P, 0j = pvj , 

, ( 6 . 2 ) 

fj = P v j v k + pSjk, j = 1, 2, 3. 

The scalar function rj = rf(9) is called an entropy function and scalar functions 
q k (9 ), k — 1, 2, 3 are called entropy flux functions, if they satisfy: 

Very (9)V e f k (9) = W e q k (9). (6.3) 

Since the the right hand side of (6.3) is a gradient of the function q k , the relevant 

integrability condition (c.f. [1], page 39) is 

( V 2 eV )(W e f k ) = (V e / fc ) T (V^). (6.4) 

If we find such an rj that is strictly convex, the change of variables 9 —» tv — 
(V eo) T will render (6.1) into the symmetric form (2.17); see [1], where 

A° = (Vfo)" 1 , (6.5) 

To see this, we apply chain rule: 

d a 9 = ('V e w)~ 1 d a w = d a w. ( 6 . 6 ) 

Substituting ( 6 . 6 ) into (6.1), we obtain 

(v^r v+£( v »/‘x v 2>r w = o. 

k=1 
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A G is positive definite if and only if rj is strictly convex. To verify the real symmetry 
of A k , we use (6.4). Multiplying both sides of (6.4) by (V^p) -1 on the left and 
right, we see A k = (A k ) T . 

We will solve (6.3) keeping the mechanical energy of classical Euler equations in 
mind. Thus, instead of 9, we will use z = (p, v 4 , v 2 , v 3 ) T as independent variables. 
We compute: 


where 


Moreover, 


V*0 


( B 1 B 2 v t \ 

y B 3 v B 2 vv T + B 4 I 3 ) ’ 


_ c 2 + p' p' _ o pc 2 +p 
1 c 2_ v 2 c 2 ’ 2 ( c 2 _ w 2 ) 2 ’ 

_ c 2 + p' pc 2 +p 

£>3 — —y -95 -04 — ~ 7 -9 • 

/oZ n 1Z pZ rt iZ 


detiy z 6) = 


(,pc 2 + p) 3 (c 4 — u 2 p') 

c 2( c 2 _ ^2)4 


> 0 , 


in fi~. We can thus compute the inverse of V z 9: 


(6.7) 


( 6 . 8 ) 


(6.9) 


with 


/ c 2 (c 2 + v 2 )E\ —2c 2 Eiv t \ 

\ -c 2 (c 2 + p')E 1 E 2 v 2p'EiE2VV T + E 2 I 3 J ’ 


( 6 . 10 ) 


E 1 


c 4 — pT 2 


Eo — 


2 2 

<r — rr 


pc 2 + p 


( 6 . 11 ) 


Based on (6.10), we will solve (6.3) using z as independent variables for convenience. 
In the z-variables, (6.3) can becomes 


V 2 p(7 fc = D z q k , k = 1,2,3, 


( 6 . 12 ) 


where 


and 


(V 2 / fc ) = 


c k = 


B 3 v k 


B 2 VkV T + B 4 e£ 


B 3 v k v +p'e k B 2 VkW + B 4 vej, + B 4 v k I 3 J ’ 


(V*0)" 1 (V,/ fc ) 

c 2 C\Vk 

~CiC 2 VkV + C 2 e k 
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Cse T k 

~C 4 ve[ + v k I 3 J ’ 


(6.13) 


(6.14) 



with 


C x = 

c — p 

C 2 = p'E 2 = 


c* — p v z 


C 3 = 

c 2 {pc 2 + p) 

/~4 — 0,2^/ 

r p'{c 2 

C4 ~ rA _ 


p'{c 2 — v 2 ) 

pc 2 + p 
_ „, 2 \ 


p'v 2 


(6.15) 


Formally, (6.12) is an over-determined system, consisting of 12 equations for 4 
unknowns. We seek solutions with the special form: 


v = v(p,y), q k = Q{p,yW , y = v 2 = v\ + vl + vj. (6.16) 


to reduce the number of equations in (6.12). Substituting this ansatz into (6.12), 
we obtain the following first order linear system: 

j f h = Qyi 

< c 2 C\p p + 262(1 — C\y)r]y — Q p ( 6 . 17 ) 

[ C 3 p P ~ ZC^yrjy = Q. 


This seems still an over-determined system. However, it is possible to derive a 
decoupled equation for Q from (6.17). We first multiply the second equation of 
(6.17) by (p 2 c 2 +p), and using (c 2 — p')C 3 — c 2 Ci(p 2 c 2 +p), we have 

(c 2 — p')C 3 r\p + 2C 2 {pc 2 + p)(l - C t y)r]y = {pc 2 +p)Q p . (6.18) 

Then, we compute (c 2 — p') x ( 6 . 17 ) 3 : 

(c 2 - p')C 3 pp - 2(c 2 - p')C A yr] y = (c 2 - p')Q. (6.19) 

We subtract (6.19) from (6.18) and substitute r\ y with Q y , using (6.15); this 
reduces (6.17) into the following decoupled system: 


V y Qy: 

2 (c 2 - y)p'Qy = {pc 2 + p)Qp - (c 2 - p')Q. 


( 6 . 20 ) 


We now proceed to solve (6.20) with the help of (6.17). First, (6.20)i gives 


V = Q{p:P) + G{p). 


Substitute this into ( 6 . 17)3 to get 


2 2 2 2 
G n = C -^-Q- C —^-Q pi 


pc 1 + p 
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( 6 . 21 ) 



or equivalently 


G p = 


1 


1 


-Q 


( 6 . 22 ) 


2 , m -^Qp, q = {c -y)Q. 

/jC z + p c z 

We observe that G depends on p only, so we have a linear first order ODE for q. 


pc 2 + p 


1 

q - -pq P = 


f(p ) 


which has the solution 

q(p) = e M {g{p) + h(y)). (6.23) 

Here (f>{p) is defined in (2.8). Substituting (6.23) into ( 6 . 20 ) 2 , and separating vari¬ 
ables, one has 

pc 2 +pdg 2 \ dh ,i ( \ fRO A\ 

Where, the first term in (6.24) is independent of y. while the second term is indepen¬ 
dent of p. Thus, m is independent of both p and y. We conclude that mn — const.. 
Thus, by integrating (6.24), we have 


q = Di(pc 2 + p) + D 2 e^ p ) \J c 2 — v 2 : 


G — - p + D 3 , 


Q = D 


pc 2 + p 
c 2 — V 2 


+ D 2 


?<t>{p) 


(6.25) 


Vc 2 — u 2 


?7 = D 


pc + p ^ e 


<Mp) 


c 2 — V 2 


+ D 2 


\J c 2 — V 2 c 2 


5-P + D3 , 


where D 1 , 1 )> and D 3 are integration constants. With K as in (2.8), one choice is 


thus, 


D 1 — c , D2 — —cK, D3 — 0, 

2 £ cKe^ p ) 

n — c p - , . 

1 H Vc 2 - v 2 


(6.26) 


The associated entropy-flux is (g , g 2 , g 3 ) T defined by 


q = 


c 2 (pc 2 +p) cKe^G) 


c 2 — V 


-Vk- 


Vc 2 — v‘ 


-rVk- 


(6.27) 


Moreover r) is strictly convex as was verified in Section 2. 

29 



References 


1. C. Dafermos, Hyperbolic conservation laws in continuum physics , Springer-Verlag, 2000. 

2. Y. Guo and S. Tahvildar-Zadeh, Formation of singularities in relativistic Fluid dynamics 
and in spherically symmetric plasma dynamics, Contemporary Mathematics 238 (1999), 
151-161. 

3. T. Kato, The Cauchy problem for quasi-linear symmetric hyperbolic systems, Arch. Ration. 
Mech. Anal. 58 (1975), 181-205. 

4. P. Lax, Development of singularity of solutions of nonlinear hyperbolic partial differential 
equations, J. Math. Phys. 5 (1964), 611-613. 

5. T. Liu, The development of singularity in the nonlinear waves for quasi-linear hyperbolic 
partial differential equations, J. Differential Equations 33 (1979), 92-111. 

6. T. Makino and S. Ukai, Local smooth solutions of the relativistic Euler equation, J. Math. 
Kyoto Univ. 35-1 (1995), 105-114. 

7. T. Makino and S. Ukai, Local smooth solutions of the relativistic Euler equation. II., Kodai 
Math. J. 18 (1995), 365-375. 

8. T. Makino, S. Ukai and S. Kawashima, sur la solutions a support compact de I’equation 
d’Euler compressible, Japan J. Appl. Math. 3 (1986), 249-257. 

9. B. Perthame, Non-existence of global solutions to Euler-Poisson equations for repulsive 
forces, Japan J. Appl. Math. 7 (1990), 363-367. 

10. A. Rendall, The initial value problem for self-gravitating fluid bodies, Mathematical Physics, 
X (Leipzig, 1991)(Berlin), Springer, Berlin (1992), 470-474. 

11. T. Sideris, Formation of singularity in three-dimensional compressible fluids, Commun. 
Math. Phys. 101 (1985), 475-485. 

12. T. Sideris, Formation of singularities of solutions to nonlinear hyperbolic equations, Arch. 
Ration. Mech. Anal. 86 (1984), 369-381. 

13. J. Smoller, Shock waves and reaction diffusion equations, Springer-Verlag, 2nd Edition, 1993. 

14. J. Smoller and B. Temple, Global solutions of the relativistic Euler equations, Comm. Math. 
Phys 156 (1993), 67-99. 

15. S. Weinberg, Gravitation and Cosmology: principles and applications of the general theory 
of relativity, Wiley, New York, 1972. 

16. Z. Xin, Blowup of smooth solutions to the compressible Navier-Stokes equation with compact 
density, Comm. Pure Appl. Math. 51 (1998), 229-240. 


E-mail address: panrh@math.gatech.edu, smoller@umich.edu 


30 



